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ABSTRACT 
In this paper the distribution of natural waves in dissipative neodnorodnonyh flat bodies. Wave 
motion is described by linear integral-differential equations. Solving this problem, we obtain the 
relationship between the velocity of the wave and its length. The problem of this kind is of great 
interest to geophysicists in the field of engineering and construction. 
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1. Introduction 


Many construction and engineering design work in dynamic conditions, is made up of 
deformable bodies having different viscoelastic properties [1,2,3]. In addition, an important role in 
the processes of wave elastic bodies due to play in signal processing tasks, especially in connection 
with the creation of mechanical resonators and filters [4,5,6]. The mechanisms by which the energy 
of elastic waves is converted into heat, are not entirely clear. Proposed various loss mechanisms [7 - 
11], but not one of them does not fully meet all the requirements. Probably the most important 
mechanisms are the internal friction in the form of sliding friction (or attachment, and then slip) 
and viscous losses in the pore fluids; The latter mechanism is the most significant in a highly 
permeable rocks. Other effects are likely, in general have less importance, are the loss of the heat 
generated in the compression phase of the wave motion by conduction, piezoelectric and 
thermoelectric effects and the energy for the formation of new surfaces (which plays a significant 
role only near the source). Therefore, development of a common methodology and algorithm for 
calculating wave fields dissipative inhomogeneous layered bodies, is an urgent task mechanics of 


deformable solids [12-13]. 
2. Formulation of the problem 


Suppose that in a Cartesian (x, y, Z) coordinate system, the beginning and the OZ axis, a sequence 


of parallel planes (fig.1) 
Z=0, Z=hi, Z=hi+hg......., Z=hithet...+hn 


The plane Z = hi + h2+ ... + hn, (n = 2) is called the n - m horizon Assume that the space between 


said planes filled with isotropic elastic medium, forming parallel layers. Layers 0 <z <h, 
characterized by a constant 4), £4)» Pp, , we call zero. Wednesday is hi + ho +... hn <z <hi + he +... hn 
+ hn +1 fills the space between the n —mand the n +1 -mhorizons, characterized by a constant 


Aes Pepi will be called the n—m layer. This will always be assumed that the adjacent layers differ 


from each other by at least one of the permanent A, LL andp. In the theoretical study of the 


processes described, we assume that within each layer of wave propagation is described by the 


usual equations of the theory of elasticity. As for the conditions at the interfaces of adjacent layers, 
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it will be assumed that when passing through them are continuous components of the vector of 


elastic displacement and stress tensor. Such contact is called hard. The paper investigates the 


dynamics of dissipative inhomogeneous planar bilayer structures. 


ell aR y = hithot...+hn 


Fig.1. Design scheme: the body in the half. 


Accounting internal friction caused by the scattering of energy in construction materials, is more 
difficult. Soft layers of multilayer structures (aggregates), as a rule, are made of materials that have 
developed rheological properties. Therefore, the dissipation of energy in the first place to be 
considered for soft layers, because it is mainly occurs during deformation of these layers. 
Mechanical systems, for which the viscoelastic properties of the elements are identical is called 
dissipative homogeneous system with different rheological characteristics - dissipative 


inhomogeneous [1,8,9]. 


The equations of motion of the deformable layer in the absence of mass forces are [1]: 


2- 


iVi+A, +f, grad did i = pi Gay239 7) 


Here U(U,,U,,,U.) -displacement vector points of the medium;g- material density; ui — moving 


parts; V ; - Poisson's ratio; 
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~ v.E. VE 
ee a J = u.= 
J is KH; 
d+v,)d—2v,) 


E— modulus operator, which have the form [9,13]: 


E,o(t)= &,,| e¢)-[ Ry (t—7)e@dz | @) 


g(t)-arbitrary function of time; R,,(t—r)-relaxation kernel; £,,-instantaneous modulus of 
elasticity; we accept the integral terms in (2) small, then the function g(t) = y(t)e”"", wherey(t)-a 
slowly varying function of time, @, -real constant. Next, using the freezing procedure [9], we note 


the relation (2) approximate species 


E,p=E,,|l-I¢ (o,)-iT 3 (op, 


wherel’,“(Wg)=|R,(c)cos@,r de, T;* (Ox) = J R,(c)sin @,tdt, respectively, the cosine and sine 
0 0 


Fourier transform of the relaxation of the core material. As an example, the viscoelastic material 
* ‘ —B;t l-a; ; 

take three parametric relaxation nucleus R ,(t)= A,e “/t- On the effect of the function 

R ,(t—7) superimposed usual requirements integrability, continuity (exceptt =T), signs - 


certainty and monotony: 
<0, 0< | R(@ar< 1. 


0 


dR(t) 


R>0 
: dt 


u- medium displacement vector j-th layer. 
On the border of the two bodies can specify two types of conditions: 


1. In the case of hard contact in the interface is put the condition of continuity of the relevant 


components of the stress tensor and displacement vector, i.e. 


dd) _ (2), dd) _ (2), 
y =O, ? Or = O01, ? 


(1) (2). .@ (2) (3a) 


Uy =u, uy =Uy. 


If the interface is no friction, the 
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<n) ae ee Oe Pee 
0, =Oy 3 Oy =O, =0; u =Uy; (2,6) 


2.On the free surface is placed a condition of freedom from stress, 


; (1) _. ) _. 
eS Og Fe Oy" (2,¢) 
where 
Ou Ou, Og. 
(j)) J Gj) J j 
o! ey ei Cc. u [Fes 4 
a) 0g. 
0) =1,0,4+2n,—29 = 4 
: ox oy 


3. Solving methods 


Now consider the solution of the differential equation (1) - (2) for a single layer. The equation of 


motion of the displacement is reduced to the following form: 


Ox? Oy’ Ox\ O Oy Or? 3) 
Oe: OG |) ee OO My Oe. ag, 
Ba ee + Panwa = 9: 

Ox Oy Oy\ Ox Oy Ot 


where (,, - material density. The solution is found in the form of: 
Un=Unly) eo > Fp = Vilyee rs HT 2: secs NA) 


whereU,(y) andV,(y) — amplitude integrated vector - function; k - wave number; C = Ce + iC; - 


complex phase velocity; and w - complex frequency. 
To clarify their physical meaning, consider two cases: 


1) k = kp C = Cpr +iC;, then the decision (4) has the form of a sine wave in x, whose 


amplitude decays over time; 
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2) k = kr +iki; C = Cr, then at each point x fluctuations established, but x decay. 


In both cases, the imaginary partkiorCicharacterized by the intensity of the dissipative processes. 


Substituting (4) into (3) we obtain: 
iz, "Ku, )+(Z, +H, )ik(iKU, +V/)+ p,k?C°U,, =0; (5) 
i,(Ut KU, )+ (2, +H, )ik(RU, +V,)+ p,K°C°U, =0. 


Thus, we have the equation (5) second-order two regions each. The problem is solved 
directly, without reducing the equation to the fourth order equation. All the arguments are for the 


layer. 


Private solutions of the system (5) is in the form 


Where rn — constant. A homogeneous algebraic system relativeAnandB,it has a nontrivial solution if 


its determinant is equal to zero 


Ce aPeC-C rier ACC 
= 0,(6) 
(C2, -C2,, jkr, Cee eae 


C? =(4,+220,)/p,. Cy =H,/P,.atn = 0 valueC7,u C’mare respectively the speed of 


compression and shear waves in an elastic medium [4]. Equation (6) may have four roots 


(7,.),, =tkf1-C7/C2; (,)4, =thf1-C?/C2 sn =O); F=1...4. 


As a result, we find four particular solutions of the form 
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U, ae (r, )i” 
Vy = Gn Blo on=O1. (7) 


Substituting values (rp) iat (7) findAni, Bniattn =(tn) i. 


Expressions for the displacement are as follows: 


U ik - ik _ 
ra e(%,, Jem se,(t Jews 
V — kq, kq, 
—kS, —kS, 
C 1 ery +4 c 1 en. 
sof 8) esse) 


U ik = ks 5 
= 4 : _ | e® +C,, a ar. 
V —kq, —ik 


Consequently, both the hard and to obtain a sliding contact set six boundary conditions 
which lead to six homogeneous equations with six unknowns C11, C12, C13, C14, C22, C24. To such a 
system of equations have a nontrivial solution, the determinant of the coefficients must be zero. 


This equation gives the dispersion equation for dissipative systems, where 
S,=G-C7/C*% 7; Gg, =A-C*/C*n), n=O, 


As an example, consider the problem of distribution of natural waves in the viscoelastic layer 


at the half. 
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Hard contact. Dispersion equation has the form 


1+ S57)e% (1457) ce 20% 2. o2e™ 2 oO. . 
aa _ Ss a 1 -Eq. on 1 z 
—2q,°" fs 2g," : {s +t) ; [= +b lee sw 00 2 4 2 
sy Sy 
(4S?) ee oes) oe «2 oe". o 2 4s )y 2/7, (8) 
1) z 1) -z, 2%, 1 i 
2g sy Bg og a {s + a [= + 1 |e ——s [s + i\y 
Sy S| “11 . 
eh 2. est rrr x te fw OM LLL 1. .-l 
gem. qe ere ae : aie on Gad A 
S Ss Ss 


where C —dimensionless wave numberC = kh,yi= “,/ orA+2U= 


= 2 (1-v)/(1-2v).As a relaxation kernel viscoelastic material will take a three-parameter kernel 


ABP nis 
R(t)= = Rizhanitsena -Koltunova [13], has a weak singularity, where A,a,f-parameters of 
t 


materials [13]. Assume the following parameters: A = 0,048; 6=0,05; a=0,l. Using the 


complex representation of the modulus of elasticity as described previously. 


The roots of the frequency equation is solved by Muller, at each iteration of the method 
applied by Muller Gaussian with the release of the main element. Thus, the solution of equation (8) 
does not require disclosure of the determinant. As an initial approximation we choose the phase 


velocity of the waves of the elastic system. For free waves in R, =0 phase velocity and wave number 
are valid values. In the calculations we take the following values: 

6 = pi/po= 0, 75; b= 107 n= 1. 

Consider two options for a dissipative system. In the first embodiment, the dissipative system is 


structurally homogeneous. 


Wave number Evaries 0 — 3. The calculation results are shown in fig.2.a. The dependence of 
the frequency and damping of the dimensionless wave number €lt was monotonous, and the 


character according to the same frequency and damping coefficients. In the second embodiment, 
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the dissipative system is structurally non-uniform: half-considered, the equation (8) and elastic 


parameters coincide with those adopted above. The calculation results are shown in fig.2.b. 
Frequency Dependence &lt is the same as for the homogeneous system: corresponding curves 


coincide up to 5%. Dependencies damping coefficients of Enon-monotonic. 
Of particular interest is the minimum value €a fixed damping coefficient: 
S=min (-w/k), k= 7,2,.....K, (9) 


hered — coefficient determining the damping properties of the system (let's call it a global damping 


coefficient). 


of ny 


Fig. 2.a. Changing the complex natural frequencies of the wave number .a) Dissipative 


homogeneous mechanical system. 
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For a homogeneous system coefficient 5 is entirely determined by the imaginary part of the 


first modulo complex frequency. For heterogeneous systems as coefficient & can act imaginary 
parts of both the first and the second frequency depending on their values. "Turn the Tables" 
occurs when the characteristic value of & the closest at this value of the real parts of the first and 


second frequencies. 6 ratio at a specified characteristic value has a pronounced maximum. 


The sliding contact. Dispersion equation is similar in form to equation (8). All parameter values 
coincide with those adopted above. Figure 3a and b shows the dependence of the frequency and 
damping coefficients of the wave number &€, respectively, for structural and homogeneous and 
heterogeneous systems. These results confirm earlier findings. Changing the parameter fromso that 
essentially dependent coefficient 5, can be achieved by varying the geometric dimensions of the 


members without altering their mechanical properties. 


- 2 3 4 
Fig.2. b Change of complex natural frequencies of the wave number. 
b) Dissipative heterogeneous mechanical system; 
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a) Homogeneous dissipative mechanical system. 


4 
b) the dissipative heterogeneous system. 
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Fig.3. Changing the complex natural frequencies of the wave number. 
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This opens a promising possibility of effective control damping characteristics of heterogeneous 


viscoelastic systems by changing their inhomogeneous systems with close frequencies. 


Analysis of Stoneley waves and Rayleigh. Consideringthehalf- 
layerandwithdifferentcoefficientsofPoisson(vi=0.3; v2=0.35; rheological properties and other 
parameters coincide with those adopted above), we see that the Rayleigh wave velocity and 
attenuation wave velocity for a layer less than Stoneley wave velocity at the boundary layer section 


- half-wave if the latter exist (fig.4). 


Fig. 4 Change C* and é the wave number. 


Results of the study of structurally homogeneous and heterogeneous systems are practically the 
same. At y/8> 1 The speed of the shear waves layer above the speed of shear waves in the half, 
then there is only one form of vibrations. The results of calculations for the hard and sliding contact 
with the original data 9 = 1.44; y = 8; vi = v = 0.25 (Other parameters coincide with those adopted 


above) are shown in fig. 5. 
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Analysis of Love waves. Consider the problem of the distribution of the Love waves in two- 


layer media, taking 


Un = [u(y),Wn (y)]e 14 * 2), (10) 


= a’w 


The wave equation (1) takes the form: = Hn ae + by? | — p,0°w, =0; n=1,2 


(11) 


-a*-10 


Fig. 5. Change Cr* and C;* the wave number ¢. 


Suppose that the boundaryy = hfree of loads. The solution of equation (9) in the form: 


paeZ 1 vd 0 


u=0, v=0; Wie Wye. (12) 
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whereW; = Asin (B1y) +Bcos (B1y), 0 <y<h 


W3=Cexp (B2y), y< 0; 
2 2 
(40) (40) 
Be ee es a: B eres a 
"GY ~  (Cp)5 


Ce/Co -C/Co1 


1.5 3,0 


Fig.6. Dispersion curves Lava stress waves in a layer of visco elastic half-space(a*= a H; C*= C/C,). 


The solution of equations (10) selected from the condition limits at infinity:Re= (B2) > 0. The 
integration constants A, B, C associated with the boundary conditions in the plane y = h and the 


conditions of continuity of stresses and displacements in the plane y = 0. 
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To satisfy the boundary conditions must be equal 


cos(Z,h) sin (ZA) O 
O 1 —1 = 0. 
Lt, B O LM, By 


In the numerical calculation will take p1/ w2 = 0,1; 8 = 0.75, other parameters coincide with 
those adopted above. The calculation results are shown in Fig. 6. The results show that when there 


Stoneley waves or Love waves, the observed effect does not occur. 


As a second example, consider the dissemination of its own waves in a flat layer being in a 


deformable (viscoelastic) medium (Figure 7). 
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Fig. 7. The dispersion curves of phase velocities Xp, , Xp and the attenuation factor V7, , 
ZX 12 for dural strip in contact with the medium acrylic 


Cp1=5400 m/c; v 1=0.35;Cp2=2300 m/c; v ,=0.35; Cs1=1311M/c; p1 =0.138 Kr/m?; 


Cs2=1311 m/c; p2 =0.126kr/m?. 


The calculation results are presented in Fig. 7. Dependence of frequency on € was the same as for a 
homogeneous dissipative system: the corresponding curves coincide up to 5%. As for the damping 
coefficients, their behavior has changed radically: the dependence w) ~ € became nonmonotonic. 


Global damping coefficient of the specified characteristic value of € has a pronounced maximum. 


4. Conclusions 


e In the course of solving the problem of wave propagation in a dissipative - inhomogeneous 
media found non-monotonic dependence of the damping rate of the physical - mechanical 


and geometric parameters of the system. In dissipative heterogeneous environments according 
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to the phase velocity and the damping rate of the geometric and physical - mechanical 


parameters of the system turned out to be non-monotonic; 


On the basis of the numerical results obtained revealed that the possibility of separation of 
thin-walled structures of the soft layer and the effect of the resonance to speed on the size of 


the site contacts. Also taking into account the material of viscous properties of 15 - 10% 


increase in the value of the phase velocity; 


Detected that higher phase waveforms and torsional expansion exceeds the highest possible 
speed (S) waves in an infinite medium, the group velocity never exceeds S. also found that the 
group velocity of a nondispersive medium 10 - 15% higher than the comparison dispersion 


medium. In other words, the shape of the pulses as they propagate not remains unchanged as 


homogeneous elastic bodies. 
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